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In the thermodynamic limit the ratio of system size to thermal de Broglie wavelength tends to 
infinity and the volume per particle of the system is constant. Our familiar Bose-Einstein statistics 
is absolutely valid in the thermodynamic limit. For finite thermodynamical system this ratio as 
well as the number of particles is much greater than 1. However, according to the experimental 
setup of Bose-Einstein condensation of harmonically trapped Bose gas of alkali atoms this ratio 
near the condensation temperature(T c ) typically is ~ 32 and at ultralow temperatures well below 
T c a large fraction of particles come down to the single particle ground state, and this ratio becomes 
comparable to 1. We justify the finite size as well as ultralow temperature correction to Bose-Einstein 
statistics. From this corrected statistics we plot condensation fraction versus temperature graph. 
This theoretical plot satisfies well with the experimental plot(A. Griesmaier et al..,Phys. Rev. Lett. 
94(2005)160401). 

PACS numbers: 05.30.Jp, 03.75.Hh, 05.40.-a 



Let us consider a noninteracting many particle system 
of Bose gas containing N alkali atoms. The system is 
in equilibrium with its surroundings at temperature(T). 
The mass of a particle is m. The length and the volume 
of the system are L and V(= L 3 ) respectively. Thermal 

de-Broglie wave length of a single particle is At = yf^P' ' 
where k is the Boltzmann constant. Average separation 
of particles is I = (j^) 1 ^ 3 - In the classical limit ^> 1, 
i.e. -^tt ^ ./V 2 / 3 . At sufficiently low temperatures when 



function e flN / kT so that Z{N) is maximized at a suitable 
negative value of fi. Thus we can approximately write 
J2 N < Z(N')e» N ' '/ kT = Z{N)e» N / kT AN', where AA^' is 
the width of the Z(N') distribution about the maximum 
Z{N). fi is called the chemical potential. In the thermo- 
dynamic limit A — * oo and length scale L — > oo. From 
the first condition of the thermodynamic limit we can 



write 



AN' 



N 



0. 



Using *f 



~ 1, the gas becomes degenerate and quantum cor- 
rection is necessary. Bose-Einstein condensation(BEC) 
occur at the onset of this degeneracy. So, at the conden- 
sation temperature k J^ ~ N 2 ^ 3 . However the condition 

of thermodynamic limit(L/AT — > oo) is also valid at this 
temperature. For finite system the thermodynamic limit 
is realized as L/Xt 3> 1. However, at sufficiently ul- 
tralow temperatures L/Xt can be comparable to 1. In 
this situation, the Bose-Einstein(B-E) statistics which is 
valid in the thermodynamic limit needs a correction to 
the statistics. 

The prescription of the derivation of B-E statistics 
is as follows [lj. The partition function for this system 

of N particles is Z(N) = J2r e -( niei+n2e2+n3e3+ ---^ kT , trapped Bose gas is L 



Using the second 
■> oo) we can write 



we can arrive at the 
relation InZ = -fiN/kT - ^ ln(l - e -^'^l kT ). The 
average no. of particles in the single particle state(i) 

is n l — g- e _ ( „ iei+ „ 2e2+ „ 3C3+ .. 0/fcT . 
condition of thermodynamic limit (L 
6i+ ,V £ ' — ► (since a ~ ^-h)- From this condition we 
can write n~i — —kT d j^ z . So, in the thermodynamic 
limit we have the B-E statistics as n; = (ei _ >j) 1 /fcr — j. 

The usual condition of thermodynamic limit is not 
properly satisfied in the case of the experimental setup of 
BEC of 3-d harmonically trapped Bose gas@, d, 0, U @j- 
In the semiclassical approximation the number density 
at the distance(r) from the center of the trap is[7fl n(r) ~ 

2 2 

e - m 2kT j where ^mtu 2 r 2 is the trap potential and ui is 
the angular trap frequency. The length scale of this 3-d 

Putting this expression 



where q is single particle energy level and n, is the 
number of particles in the «th single particle state. For 
Bose gas rij = 0, 1, 2, 3... R represent all possible number 
distribution({rii}) subject to the constraint ni ~ 

N. For this constraint the evaluation of Z{N) without 
approximation is impossible. Since Z(N) is a rapidly in- 
creasing function we choose a exponentially decreasing 
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of length in k J^ ^> A 2 / 3 we get the condition of clas- 
sical limit as W- N 1 / 3 , condition of thermodynamic 
limit(L/AT — > oo) as W- — > oo and we get the condensa- 
tion temperature for this system as T c ~ i^A 1 / 3 . More 
precisely the thermodynamic limit for this system is re- 
alized as uj — > 0,N — > oo and Nu> 3 — constant. Ex- 
act calculation with thermodynamic limit show thatQ 

l - v / HS r and B 9 Tc .= ^& /3 - With 

the consideration of thermodynamic limit, total number 



2 



excited particles at T < T c are N e = [|^] 3 C(3) and the 
number of condensed particles are N = N[l — (y-) 3 ]. In 
the experimental setup the typical value ofQ N, uo and 
m are of the order of 50000, 2645s -1 and 52 amu. With 
these experimental parameters the length scale of the sys- 

56 x 10 -4 mm. and 



tern at T c is L 



/ 2C(4)fcT 
muj y £(3)hu 

\t c = 2.88 x 10~ 4 mm. and their ratio is L/Xt = 32.2. 
However, at sufficiently ultralow temperatures (T ~ ^) 
well below the condensation temperature, a large fraction 
of particle come down to the ground state and the length 

scale of the system becomes ~ \ ~ 6.76 x 10 _4 mm. 
Here we see that at ultralow low temperatures well below 
the condensation temperature, the thermal de Broglie 
wavelength becomes comparable to the system size. At 
these ultralow temperatures(T ~ %) the usual theory 
of statistical mechanics of finite system is not (properly) 
valid. So at ^ < T < T c we seek a ultralow temperature 
as well as finite size correction to B-E statistics. 



To quantify the correction arising from ultralow tem- 
peratures, let us start from Tsallis type of generalized 
Bose-Einstein statistics as n, = t , 



ttddq- 

where q is a hidden variable [10|, [Tl|. As q — ► 1 
get back B-E statistics. The relative probability that 
E be the total energy of a system is given by Boltz- 
mann factor e~ E ^ kT . In Tsallis statistics this factor 
is replaced by tttt— T^fpfmw^) ■ As 9^ 1> we S et 



-1 

we 



(l+(q~l)E/kT) 1 /(l~^ ■ ^ q 

back Boltzmann statistics. However, Tsallis statistics 12 
( n + (g-i)ff/fc T )i/(3-i) ) is a PP lied to equilibrium [UH!, Hi, 
[l5| as well as nonequilibriumflfl [l?], EH systems. Tsal- 
lis statistics is rederived as one of the superstatistics [l9j 
of a nonequilibrium system. In this theory the hidden 
variable (q) of Tsallis statistics is equated with system 
parameter and q is no longer a variable. In the theory of 
superstatistics[19] 2/(q — l) is redefined as effective no. of 
degrees of freedom. In the theory of dynamical founda- 
tion of nonextensive statistical mechanics, this effective 
no. of degrees of freedom is equated as[16( (3 — ?)/(<?— 1). 
However, for finite equilibrium system we equate q with 
a system parameter. We equate q — 1 with At / ' L so that 
at the thermodynamic limit(L/AT ^> 1) of finite system 
we can go to the usual Bose-Einstein statistics. Fig.l 
shows experimental and theoretical plot of condensation 
fraction with temperature. In this figure we see that a 
negative shift of condensation fraction is necessary to sat- 
isfy the theoretical plot with the experimental data. In 
the theoretical plot there is finite size correction and the 
correction of two body interaction. We shall see that the 
generalized B-E statistics with (q — 1) cx |^ will give 
rise to a significant negative shift of condensation frac- 
tion with temperature. This significant shift along with 
finite size correction and the correction of two body in- 
teraction might satisfy the experimental results. So for 
the trapped Bose gas we shall start from the general- 



ized statistics m 



where a is an 



arbitrary constant. We shall determine this a from the 
experimental result. 

For a system of Bose gas in isotropic harmonic trap 
single particle energy levels are ej = (I + j)hoj, (j = 
0, 1, 2, 3, ...). If T c be the condensate temperature then at 
T < T c the chemical potential fi = ^Huj. So at T < T c the 
number of particles from the generalized B-E statistics 
follows as 



(1) 



J ]at 
at 2 J 



where t — Obviously, in the thermodynamic limit 

of this finite system t ^> 1. Since for t ~ 1 the thermal 
de Broglie become comparable to the system size this 
expression of nj is not valid for t < 1. For a 3— d isotropic 

harmonic oscillator the density of states is - + ^ 3+2 ^ . So 

total no. of particles in the excited states at t < t c = ^j- 
will be 



N, = 



U 2 + 3j + 2) 



1 



3=1 



[1 



J ]at 
at 2 ' 



1 



(2) 



In this summation the volume term of the density of 
states(i.e. the term j ' /2) will dominate over the sur- 
face term (i.e. the term 3j/2). This surface term gives 
finite size correction. The third term contribute insignif- 
icantly to the calculation of number of excited particles. 
Converting the summation into integration we have 



f + 3j 



1 



o 



[1 



J ]at 



1 



dj. 



(3) 



For t » 1 the term [1 + -^i] — 1 in the denominator of 
the above equation can be approximated as [1 + ~ 

At t 3> 1, considering the significant correction terms 
from the above equation(3) we can write 



N e = 



°° J 2 



1 



2 e^V* - 1 



dj 



3j 



1 



2 &»'/* - 1 



dj 



2at 3 



t 3 C(3) 



°° j 2 



J/t 



2 {eJ/t - I)- 



dj 



3* 



C(2) 



(it 2 



C(4) 



(4) 



Till now we did not consider the particle particle inter- 
action. To achieve BEC it is necessary to take a very 
dilute gas. The gas being very dilute there should be 
a correction term in the expression of N e due to two 
body scattering. This correction within Hartree-Fock(H- 
F) approximation has been discussed in According 
to H-F approximation, the correction term to the above 



N e is 3 x 1.326 



y/h/r. 



7/6 



4.93 



=i 7 / 2 , where a is 



3 



the s-wave scattering length. So more corrected expres- 
sion of number of excited particle at 1 < t < t c would 
be 



V 2 

N e = < 3 C(3) + — C(2) + 4.93 



t 7/2 + — C(4) (5) 
h a 



At T = T c all the particles will be in the excited states [7], 
HHH. So at T = T c or at t c = ^ the number of excited 
particle will be equal to the total number of particles. 
So, 



fi/ 2 

* c 7/2 + ^C(4) (6) 
h a 



V 2 

JV = t?C(3) + -= £ C(2)+4.93 



In the thermodynamic limit [3, [g| iV e = t 3 ^(3) and the 
condensation temperature T is such that@, [§| t a = 



N 



1/3 



where t Q = ^ 



Comparing this expression 

of £ Q and t c of equation (6) we see that t c < t Q and there 
is a shift 5t c = t c — t Q of condensation temperature due 
to the inclusion of finite size correction term, correction 
term due to two body scattering and due to the ultralow 
temperature correction of B-E statistics. At t < t c the 
fraction of number of particles in the ground state is 
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FIG. 1: Condensation fraction (jf) versus temperature 
(-p)plot. The thick line follows from equation(7). The dashed 
line corresponds to the thermodynamic limit and excludes all 
the correction terms in equation(7). The dotted line corre- 
sponds to the finite size correction and the correction due 
to interaction and exclude the ultralow temperature correc- 
tion term in equation (7). All the theoretical curves are 
drawn according to the following experimental parameters. 
The dotted points are experimental points of Bose-Einstein 
condensation of 52 Cr where [j| N = 50000,^ = 2645s -1 , 
m — 52a.m. u., a — 105as and T ~ 700nK 



No ^ N-N, 
N 



1 



t 



N ' L t 

If 2 n (\i 2 

[£_ C(2 ) + 4 .932^=t 7 / 2 + 2L C (4)]/[^C(3)] 



(7) 



However in the thermodynamic limit the expression of 
this fraction would be@,ll| [^]t-l = 1 - [£] 3 - Now, 
from equation(7), due to the correction terms we get the 
fractional change in condensation temperature as 



5T t 



St c 



■it?. 



= X = [[[[^-[^rC(2)+ 4.932 



T, 



to 



tV 2 



Gt 



£ C(4)]]/C(3)] 1/3 - [A/C(3)] 1 / 3 ]]/[A/C(3)] 1 / 3 



1 3f 2 

-_[_£ C (2)+ 4.932 



— C(4)C(3)]/iV 
a 



tl' 2 



(8) 



Putting t c ~ [^(l)] 1 ^ 3 m the above equation we get 



He 

77 



C(2) 



2[C(3)]2/3 

2C(4) 
a[C(3)]V3 



A~ 1/3 - 1.326 



.JV-V3 



-.N 



1/6 



(9) 



From the first term of the above equation(9) we get 
the T r shift due to the finite size correction as 1221 



st f ~ 3 



C(2) 

t„ - 2[C(3)] 2 /^ 

particles we get 



N- 1 / 3 = -.728A- 1 / 3 and for 50000 



st f ~ 



-1.97%. From the second 



term of the above equation(9) we get the T c shift due 
to the correction of two body interaction as @, |2C| 



-1.326- 



a = 105a B and for N=50000 [9] 

k(3) 



6.61% for uj = 2645s" 



For this setup t = 
1 / 3 = 34.65. From the third term of the above equa- 
tion^) we get the T c shift due to ultralow temperature 

m H 2C(4) jy-1/3 = 



correction of B-E statistics as 



*[C(3)] 2 



-5.19% for 50000 particles. According to the experi- 



ment ^ should be 10% and since ^f L 



stf- 



+ 



St' 



St 



to to 

3.6549 However, since these shifts are 



, we get a 

not negligible, in equation (8) we have to take higher 
order terms into account. As the corrections are not neg- 
ligible we should not put t c = t Q in equation (9). Since 
t c < t a , the percentage t[~ s and t\ nt shift we calculated 
will be lowered. Using equation (7), from the numeri- 



cal plot of with temperature we get 



= -1.93%, 



st 



l = —5.75%. From equation (7), in the numerical 
plot 10% t c shift will be achieved if we put a — 1.48. 
Since in the experimental setup, below the condensa- 
tion temperature ^fr ~ t < 34, we introduced corrected 
statistics (rij — ^ + j 1 ^ t - ) apart from B-E statistics 



4 



( n j = e j/'-i ) which is absolutely valid when t — > oo. 
Well below T C) as t approaches to 1 the B-E statistics 
needs more corrections. Below T c as a comparative study 
of the two statistics, let us calculate the ratio of specific 
heat from the two statistics. As a comparative study we 
can take ideal gas and we disregard the surface term in 
equation (3) to get 



T(ait — 3) 
r(mt) 



(10) 



As t > 1,0, [I] N e (t) -> i 3 C(3) as expected from B-E 
statistics. The total energy of the system would be 



E(t) ~ (fiw) 



o 2 [l + ^-l 



dj 



i=i 



' T(ait) 



(11) 



Now the specific heat would be C v (t) 



k d 
hu dt 



2j 



alog[l 



E(t) 



7TP 



dj. 



(12) 

As i > 1, E(t) — > 3haj((4)t 4: as expected from the B- 
E statistics 0, |1] . Let us denote the ratio of specific heat 
calculated from the corrected B-E statistics and from the 
B-E statistics as 



r\t) 



C y 



fcl2C(4)i 3 



(13) 



Obviously as t — > oo, r' — ► 1. With a = 1.48 numerical 
plot of r(t) at very low temperatures (1 S, t < t c ) is shown 
in the FIG. 2. In FIG. 2 we see that at high temperature 
t 1, the specific heat behaves well according to our 
familiar T 3 law. At to = 34.65, the numerical value of r 1 
is 1.15 and at t — 15, the numerical value of r 1 is 1.39. So 
at t — 15 there is 39% error in the specific heat. In this 
figure we see that this specific heat law dose not holds 
fairly well at t < 6. From this figure we also see that to 
achieve the thermodynamical limit, uj — > is no longer a 
necessary criterion. If u ~ j^- , the system (3-d isotropic 
harmonically trapped Bose gas) will correspond to ther- 
modynamic limit with ~ 39% error in the specific heat. 
Similarly we can say that to achieve the thermodynamic 
limit of a system of particles in a box the volume need 
not necessarily tend to infinity. In this way we can also 
estimate the minimum volume of a thermodynamical sys- 
tem. In the FIG. 2 we see that specific heat for this finite 
system becomes negative at t < 5. At this range of tem- 
peratures L < Ay and the theory of statistical mechanics 
is not valid. In spite of that, due to nonextensivity the 
appearance of negative specific at these temperatures is 
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FIG. 2: (r'(t) 



and Temperature (t 



) ,where 



C v follows from equation(12).The doted line shows the theo- 
retical value of the ratio r'(t) in thermodynamic limit when 

LU -> 0. 
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FIG. 3: Condensation fraction {-#) versus temperature 
(-p)plot. The thick line follows from equation(7). The dashed 
line corresponds to the thermodynamic limit and excludes all 
the correction terms in equation(7). The dotted line corre- 
sponds to the finite size correction and the correction due 
to interaction and exclude the ultralow temperature correc- 
tion term in equation (7). All the theoretical curves are 
drawn according to the following experimental parameters. 
The dotted points are experimental points of Bose-Einstein 
condensation of 87 Kb where[li| N = 40000,lj = 1140s" 1 , 
m = 87a.m.u.,a = 90a B and T a = 280nK 



not surprising theoretically [23j| and experimentally [13]. 

For smaller uj the system size(L = ^ y ^^^J ) is 

larger and ^jP- is smaller. In this case the the ultralow 
temperature correction would be smaller. In the experi- 
ment [25[ ui = 1140s" 1 smaller than that of the previous 
case. For this experimental setup the numerical value 



of 



<5/ 



= -2.1% and 



to 



= —5.5%. From the experi- 



5 



mental data as shown in FIG 3, we see that ^— ~ 8%. So 

= —.4% and it corresponds to a = 5.5. Due to the 
larger system size the ultralow temperature correction is 
smaller as we see in the FIG 3. 

In the thermodynamic limit of a finite system L/Xt S> 
1. Bose-Einstein statistics is absolutely valid at this limit. 
However for trapped Bose gas at T < T c this ratio is[9( 
~ 32 and at temperatures well below T c this ratio is com- 
parable to 1. For this reason we expect a correction to 
B-E statistics. The correction to the statistics depends 
on this ratio. As a correction we introduce a nonexten- 
sive type of statistics (n.; = * akT — ) . Smaller 



the ratio larger is the correction and it is evident in FIG 1 
and in FIG 3 of the BEC experiments. Although the ex- 
perimental uncertainties are larger with respect to these 
theoretical corrections yet the ultralow temperature cor- 
rection might be of our interest. However, this ultralow 
temperature correction is valid only for finite system. 
How a is to be determined theoretically remains an open 
question. 
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